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It is well known that the dynamics of a one-dimensional dissipative system driven by the Ginzburg- 
Landau free energy may be described in terms of interacting kinks: two neighbouring kinks at 
distance i feel an attractive force F{t) « exp(—^). This result is typical of a bistable system 
whose inhomogeneities have an energy cost due to surface tension, but for some physical systems 
bending rigidity rather than surface tension plays a leading role. We show that a kink dynamics is 
still applicable, but the force F{L) is now oscillating, therefore producing configurations which are 
locally stable. We also propose a new derivation of kink dynamics, which applies to a generalized 
Ginzburg-Landau free energy with an arbitrary combination of surface tension, bending energy, 
and higher-order terms. Our derivation is not based on a specific multikink approximation and the 
resulting kink dynamics reproduces correctly the full dynamics of the original model. This allows 
to use our derivation with confidence in place of the continuum dynamics, reducing simulation time 
by orders of magnitude. 

PACS numbers: 05.45.-a, 05.70.Ln, 02.30.Jr 
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I. INTRODUCTION 

The continuum description of a physical system requires to define a suitable, coarse grained order parameter /i(x, t) 
and to build a free energy T if the system is at equilibrium, or a partial differential equation (PDE) obeyed by h if 
the system is out of equilibrium. In some cases, the PDE itself can be derived by some free energy. This is surely the 
case for a system relaxing towards equilibrium (think to a phase separation process Q) , but it may also be true for 
pattern forming systems, in which case is a pseudo free energy (also called Lyapunov functional 0). 

Typically, is made up of a potential part U{h), which is the energy density for an homogeneous state, plus 
a part which accounts for the energy cost of the inhomogeneities of the order parameter. The simplest way to 
weight spatial variations of /i(x, f) is to consider a term proportional to (V/i)^. This surface tension term appears 
in completely different contexts, from magnetism to surface physics. In the former case, the misalignment of spins 
produces an energy cost which is proportional to the gradient square of the magnetization Q. In the latter case, 
if the energy of a surface of local height /i(x) is proportional to the total extension of the surface, we simply get 
S = f dxi/l + iyhy ~ 5*0 + ^ / dx{Vh)‘^, where Sg = f dx is the area of the system. 

If surface tension combines with a double well potential U{h), which accounts for the existence of two macroscopic 
stable states, F = J'gl is called Ginzburg-Landau free energy and it plays a relevant role in the theory of phase 
transitions and phase ordering. In one dimension, a simple description of energetics and dynamics can be given in 
terms of kinks (^. A kink hy-{x) is the simplest non-homogeneous state which interpolates between the two minima of 
the potential, and it has two main features: it is a monotonous function, and it is localized, i.e. its derivative is 

exponentially small except in a finite size region. The explicit expression of a kink for a specific potential, see Eq. 0, 
h^{x) = ta.nh{h^x/, make both properties obvious. 

The reason why kinks play a major role derives from the possibility to describe h{x, t) as a sequence of kinks and, 
finally, by the possibility to describe the continuum dynamics in terms of an effective dynamics between kinks, which 
act as fictitious, interacting particles. In poor terms, kinks have an attractive force which decreases exponentially 
with their distance: the attractive force implies instability and coarsening; the exponential dependence with distance 
implies coarsening is logarithmically slow. 

In spite of the widespread importance of J^gl, we should not come to the wrong conclusion that its form is universal. 
This caveat is particularly appropriate if bending rigidity is important: soft matter and biophysics, dealing with 
membranes @ and filaments are a relevant example. This fact, the relevance of bending rigidity with respect 
to surface tension, is not purely phenomenological. On the contrary, it has been recently derived rigorosuly from 
an hydrodynamic model [7|. According to this model, the energy cost of inhomogeneities is proportional (in a one¬ 
dimensional model) to the squared second spatial derivative of h, {hxxY, rather than to the squared derivative, {hxY- 
This modification is of paramount importance, because kinks are no longer monotonous functions and this fact will 
be seen to change drastically their dynamics, which turns out to be frozen. 

The goal of our manuscript is twofold: first, we extend J^gl to a free energy which depends on surface tension, 
bending and possibly higher order terms. Second, we reconsider the problem to pass from a continuos formulation of the 
dynamics to a discrete description in terms of kinks, proposing a new approach. A detailed numerical comparison with 
continuum dynamics reveals that standard approaches where the order parameter profile is written as a superposition 
of kinks fail to reproduce quantitatively exact dynamics. Instead, our new approach is quantitative. 

The paper is organized as follows. In Section II we define the various continuos models and in Section III we 
give a simple derivation of known results. In Section IV we propose a new derivation of kink dynamics and compare 
numerically different approaches. In Section V we discuss the stability of steady states and in Section VI we summarize 
the results. 


II. CONTINUOUS MODELS 


As explained in the Introduction, a good starting point to introduce the dissipative dynamics of interest for us here 
is the Ginzburg-Landau free energy. For a scalar order parameter in a one-dimensional system, 

J-G. = J dx (^^hl + U{h)^ , (1) 

where U{h) is an arbitrary symmetric double well potential, with two equivalent minima for h = ±/i„, which are the 
ground states of the full free energy. If U{h) = UgU(h), rescaling space we obtain 


F GL — \/FlUg I dx 


= GqF gl- 


( 2 ) 
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In the following the energy scale eg will be set equal to one, while we don’t rescale to one for pedagogical reasons. 
Furthermore, for definiteness, in this Section we consider a standard quartic potential, 17(h) = —h^/i^/2 + h‘^/4. 

The free energy J^gl is the starting point to study the dissipative dynamics when the system is relaxing towards 
equilibrium. When studying dynamics the existence of conservation laws is of primary importance and two main 
universality classes exist, depending on whether the order parameter, h(x, t), is conserved or not. In the two cases we 
obtain, respectively, the Cahn-Hilliard (CH) and the Time Dependent Ginzburg Landau (TDGL) equation, 

r -p 

dth(x,t)= - ^ = h,, + hlh-h^ TDGL, (3) 

da 

dth(x,t) = = -dxx(hxx + h^h - h^) CH . (4) 

oh 

In both cases, it is straightforward to show that 


dJ~ gl 
dt 


f , ST dh „ 


( 5 ) 


Equation © typically describes phase separation in a magnet, because in this case relaxation dynamics does 
not conserve magnetization. Equation (jT]) can instead describe phase separation in a binary alloy, where matter is 
conserved. Here we will focus to a so-called symmetric quench, where the average value of the order parameter is 
zero. 

In the above two cases, TDGL and CH equations, the overall picture of dynamics is well known. The solution h = 0, 
corresponding to the disordered or homogeneous phase, is linearly unstable, as easily seen by a stability analysis. In 
fact, if h(x,t) = , to first order in e we find 


cr(g) = 


hl-q^ 
hW - 


TDGL 
CH ’ 


( 6 ) 


so that the homogeneous solution is linearly unstable for small q. Because of such instability, small regions of the two 
phases h = appear, separated by kinks. A kink is a steady solution of TDGL/CH equations which connects the 
two minima of the potential U(h) for x —>■ ±oo. For the standard quartic potential, such solution has the simple form 


h{x) = ±h^{x) = ±/i,j,tanh . (7) 

More generally, TDGL/CH equations have periodic solutions of arbitrarily large wavelength which can be thought 
of as superpositions of kinks {h]^(x)) and antikinks {—h^(x)). These kinks feel an attractive interaction, and when a 
kink and an antikink meet they annihilate, therefore leading to an increasing average distance between the remaining 
ones (coarsening process). In an infinite system this process lasts forever, but in one dimension it is logarithmically 
slow. 

The above picture is well known and goes back to works by Langer Q and Kawasaki and Ohta 0 . The main idea 
is to write h(x, t) as a suitable superposition of positive and negative kinks, getting a set of discrete equations for 
their positions Xn{t). This approach will be discussed in the next Section. First, we need to show how this picture 
should be modihed if the surface tensione term (K^) in the GL free energy is replaced by a bending term (h^^). 

If bending rigidity dominates over surface tension, the Ginzburg-Landau free energy should be written 


T GL4 — / dx 


U(h) 


and Eqs. dSEl) are modified as follows, 

dth(x, t) = -hxxxx + hlji - TDGL4, 
dth(x, t) = -dxx {-hxxxx + hl^h - h^) CH4, 


( 8 ) 

(9) 

( 10 ) 


where the label ‘4’ highlights the replacement of a second spatial derivative with a forth spatial derivative. In its turn, 
the linear spectra ([6]) should be replaced by 


J hl-q^ TDGL4 

1 hW - 9® CH4 ’ 


( 11 ) 


showing that the homogeneous state is still unstable for large wavelength fluctuations. 
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FIG. 1. Plot of kinks appearing in TDGL/CH (dashed line) and in TDGL4/CH4 (full line). In the latter case, the tail 
continues to oscillate around ±/im, but its exponential decay allows to make visible only the hrst two oscillations. 


In spite of these similarities, the study of steady states is not straightforward as for TDGL/CH, where it essentialy 
boils down to solve the problem of a particle of coordinate h in the potential V{h) = —U{h). Steady states are now 
determined by the time independent equation 

- Kxxx - U'{h) = 0. (12) 

The forth order derivative introduces new classes of kinks, because fixing the conditions h{x —>■ ±oo) = is no 
longer sufficient to uniquely determine a solution. According to Ref. kinks can be labeled by their number of zeros, 
i.e. the number of points where the kink profile vanishes (Eq. ([7]) shows that for TDGL/CH kinks this number is equal 
to one). The asymptotic behavior, i.e. the limiting behavior of hj^{x) for large |x|, is determined by the linearization 
of Eq. (fl^ around h = h^, 

hxxxx = -U'\h^){h - h^). (13) 

It is easily found that h{x) = + R{x), where the tail R{x) is given by 

R{x) = Acos{kx + a)exp{—Kx), (14) 

where k = {U"/\/2, while the amplitude A and the phase a are undetermined within the linear theory. The 
exact shape of kinks for TDGL and TDGL4 models is plotted in Fig. [U where we limit for TDGL4 to the kink with 
only one zero. 

A similar picture, oscillating kinks and kinks with more zeros, emerges in other PDFs, e.g. the convective Cahn- 
Hilliard equation |10l |. In both cases there is no evidence of such multihump kinks during dynamics, which lead us to 
assume they are dynamically irrelevant. Therefore in the next Section we are studying kink dynamics assuming kinks 
which cross the horizontal axis only once. 


III. KINK DYNAMICS MADE SIMPLE 


The following, semiquantitative treatment of a profile simply consisting of the superposition of a negative and a 
positive kink allows to grasp the relation between the kink tail R{x) and the kink interaction. In order to get a result 
as general as possible, we consider an energy functional which is the sum of a symmetric double well potential (as 
before) plus arbitrary quadratic terms, whose only constraint is to satisfy the symmetry x —> —x. Its most general 
form is 


T = 


dx 


U{h) 


M 


(15) 


where a 2 i are constants and the notation {dl.h) means the I— th order spatial derivative of h. We have also introduced 
the factor ^(—1)® so as to get rid of it when evaluating the functional derivative, according to the relations 


M 

=U'{h)-Y,a2,dl^h, 

=U'{h)-C[h]. 


5 ^ 

Sh 


(16) 
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The model we are going to analyze is a nonconserved, purely dissipative model, where dynamics is driven by T 
according to the relation dth = —(bTjbK), i.e., 


dth = C[h] - U'{h). 

If h]^{x) is the kink profile centred at x = 0, the two-kinks approximation amounts to writing 

h{x, t) = + xo{t)) - - xo(t)) - 


(17) 


(18) 


where the kinks are centred in ±Xo(t) and the constant term must be added in order to get the correct values in the 
different regions (for an kinks approximation, the constant term is more complicated, see Eq. (2.8) of Ref. @ and 
Eq. (l30l) here below). Using Eq. (fT8l) it is easy to evaluate dth, 


and its spatial integration. 


dth{x,t) = xo{hl{x + xo) + hl{x - xq)), 


/ + 00 

dxdth{x,t) = Ah^xo- 

-OO 


As for the RHS of Eq. (flTl) , while we simply have 

C\h] = C[K{x + Xo)] - C[K{x - Xo)], 


(19) 


( 20 ) 


( 21 ) 


the evaluation of U'{h) is a bit more involved. As soon as jxj ^ a, a being the size of the core of the kink, 
h^{x) ~ + ^(ja^l)], for x ^ 0 respectively. Therefore, we can approximate Eq. as follows 


h( +'1 ^ / ^k(a:-I-Xo)-I-i?(-x-I-Xo) for x < 0, 

~ \ —hj^{x — Xo) -I- R{x -I- Xo) for x > 0, 


and write, in the two cases, 


rr/fh\ ^ f U'{hy,{x + Xo)) -k U"{K{x + Xq))R{-x + Xo) 
^ ( -U'(K(x - Xo)) -k U”{K{x - xo))R{x -k Xo) 


for X < 0 
for X > 0 


( 22 ) 


(23) 


so that 


/ +00 p+ca 

dxU'{h)= / dx [t/'(/ik(x-k Xo)) — C/'(/ik(x — Xo))] 

-OO j —OO 

+ / dx [t/'(/ik(x - Xo)) -k t7"(dk(x-k xo))i?(-x-k Xo)] 

J —OO 

^+oo 

-k / dx [— C/'(dk(x -k Xo)) -k t/"(dk(x — xo))i?(x -k Xo)]. 
do 


(24) 


In the previous expression, a simple change of variable in the second line integral, x —> —x, shows it is equal to the 
third line integral. 

We can now match the spatial integration of the two sides of Eq. (El). Using Eqs. (I2()I21I24L we obtain 


^ + oo 


4d„,xo = 


dx{C[h] - U'{h)) 

dx[c[h^{x -k Xo)] - U'{h^{x -k xo))^ 
dx[c[h^{x - Xo)] - U'{h^{x - xo))) 

-k2 J dx(jj'{h^{x + Xo)) — U"{h^{x — xo))R{x + xq)^ ■ 


' —OO 
/‘ + 00 


' —OO 

/* + oo 


' —OO 

^+oo 


( 25 ) 
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Since the integrands in the second and third line vanish, we finally get 


dx(^U'{h^ + R{x + Xq)) — U"{h^{x — Xq))R{x + (26) 

= ^ - U''{K{x - xo)))r(x + xo) + O(R^). 

The quantity within large brackets in the final integral is exponentially small when |x —2:o| ^ a, so we can approximate 
the integral as the integrand value for x = xq times the extension over which the function in square brackets is non 
vanishing, i.e. a. Finally, we can write 



J dx(u' {h^{x + Xq)) — U”{h^)x — Xq))R{x + Xg)^ 


xo^^[U"{h^)-U"mR{e), 


(27) 


with £ = 2x0- In conclusion, the speed of the right kink is barely proportional to R{£), where £ is its distance from 
the left kink (the quantity in square brackets being positive, since U"{h^) > 0 and 17"(0) < 0). 

This result means that a kink exerts a force on its right neighbour at distance £, force which is proportional to i?(f'), 
where R{x) is the difference between the kink profile and its limiting value for large, positive x, R{x) = h^{x) — 

For the standard TDGL equation, h^{x) = h„,tanh(^x) and R{x) = R 2 {x), with 

R 2 {x) = —2/i„ exp(—-^x), (28) 

V2 

while for TDGL4, i?(x) = R^ix) = A cos {kx + a) exp (— kx), see Eq. (fTdll . 

We can assume that Eq. (I?fl) may generalize to any sequence of kinks located in x„(t) (with x„+i > x„), 

in = -rjT^ {U"{hJ) - U"{0)) [i?(x„ - x„_i) - i?(x„+i - x„)], (29) 


where the size a of the kink core has been evaluated as a = 2h^/h[^{0). Above equation should be supplemented 
by the constraint that two neighbouring kinks annihilate when they overlap (see details on numerical schemes in 
Appendix [B|) . 

As a matter of fact, such kink dynamics can be derived using a superposition of N kinks. 


h{x, t) = (-l)”/ik(a; - Xn{t)) 

+ ^(-l)''[/ik(x - Xfc(t)) - hj\ 

k<n (30) 

+ '^{-l)''[K{x - Xk{t)) + h^]. 
k>n 

This approach was initially used by Kawasaki and Ohta Q to study TDGL and GH equations. In the next Section 
we are going to propose a novel approach and to compare both with numerical integration of the full continuum 
equations. 


IV. IMPROVED KINK DYNAMICS 

We now provide a more general approach to kink dynamics: we don’t assume explicitely a specific “multikink” 
approximation, as, e.g., Eq. (l30l) . and we consider the general energy functional given in Eq. dig. We don’t claim 
our approach is rigorously founded: its validity (and usefulness) are rather supported by the final comparison with 
numerics. 


A. Nonconserved case 

The nonconserved case corresponds to the dynamics 

^th = -^ = J2<^2^^l^h-U\h). (31) 

i 
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FIG. 2. (Color online) Schematic of studied system with relevant notations. 



FIG. 3. (Color online) Exact dynamics and analytical approximations of the motion of two kinks for the TDGL model. Black 
squares: exact dynamics (integration of Eq. (|3}). Red full line: our model, Eq. (I39II . and Ei and Ohta’s model. Blue dashed 
line: Kawasaki and Ohta’s model. 


In Fig. [2] we show the schematic of the system. It has been drawn for TDGL4/CII4 kinks, but notations are 
generally valid. More precisely, means the position of n—th kink and the points halfway between kinks n 

and (n ± I). For ease of notation, is replaced by n± i in integrals’ extrema and h{x^^i) is replaced by 

We assume that apart from the annihilation process, which occurrs when = Xn+i — Xn ^ a, kinks retain their 
profile when moving. So, for x around Xn the previous equation can be rewritten as 

-Xn^^h = Y,a2^^fh-U'{h). (32) 



FIG. 4. (Color online) Exact dynamics and analytical approximations of the motion of four kinks for the TDGL4 model. 
Black squares: exact dynamics (integration of Eq. (O). Red full line: our model, Eq. (1431) . Blue dashed line: Kawasaki and 
Ohta’s model. 
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FIG. 5. (Color online) 

Exact dynamics and analytical approximations of the motion of four kinks for the CH model. Black squares: exact dynamics 
(integration of Eq. Q). Red full line: our model, Eq. (1481) . Blue dotted line: our model, Eq. (1501) . Green dashed line: 

Kawasaki and Ohta’s model. 



X 


FIG. 6. (Color online) Exact dynamics and analytical approximations of the motion of four kinks for the CH4 model. Black 
squares: exact dynamics (integration of Eq. (1101) 1. Red full line: our model, Eq. (1491) . Blue dotted line: our model, Eq. (1511) . 
Green dashed line: Kawasaki and Ohta’s model. 


We then multiply both terms by dxh and integrate between and : 


(•n+i 


da; (dxh)'^ = ^ a 2 i 


rn-\-4 


da; dxhdx'h -U{h^+i) + U{h„_i). 


Direct integration and integration by parts give 


Xn = 


In-l da; (dxh)^ 


E 


02 * 






k< 7 




(33) 


We stress that above result derives from one single assumption, dth ~ —Xndxh for x close to a;„. Equation (1^ can be 
further elaborated because in the region halfway between a;„ and a;n,+i we can expand h{x,t) around the asymptotic 
values ±/i„, 


h{x, t) ~ ± [h^ + R{x — Xn) + R{Xn+l — 2;)] 


(34) 
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where +/— applies for a positive/negative n—th kink. Using this notation, we finally get 

2 / « \ 2\ 


1 


Xn =- 


/_r da: 

k< 


E 


a2i 




-‘‘ii (y) (y) -(%^) (' 


where, at denominator of Eq. (I33p . we made the approximation 

^n+i 


dx [dxhf ~ 


ra+i R+oo 

dx {dxKY ~ / 


+ 217 "(/i„,) i?' 


dx {dxKf 


-i?" 



(36) 


i.e., we have assumed that close to x„ the kink profile is similar to the static profile h^x) and we have extended the 
extrema of the integral to ±oo, because is concentrated around Xn- 

Therefore, in the general case of an equation with several terms a 2 i Y 0 fhe expression of the speed of a kink is 
fairly complicated. One remark is in order: the different contributions to the RHS of Eq. (1351) are not proportional to 
R{in) and as appearing in the simple approach given in the previous Section. This point is better clarified 

by focusing on two explicit cases. 

• For TDGL, the only nonvanishing term in the summation is 02 = 1, so Eq. (1551) strongly simplifies to 


Xn = 


2 C/"(/r„ 




-R^ 


/_oo da: 

which must be compared with Eq. dMl), rewritten here for convenience: 

1 


[NEW approach] 


(37) 


Xn = 


h'M 


{U"{h^) - U"{0)) [R{in-i) - i?(^n)], [KO approach] 


where KO stands for Kawasaki and Otha [ 8 |. 

In the specific TDGL case R{£) is a simple exponential, so that 

R{£) = constant x 7?^ ( y ) • 


(38) 


In conclusion, the new approach (1571) and the old approach (1751) differ for the prefactor only. Let us work out the two 
prefactors for the explicit expression U{h) = —^ + ^. Using Eq. ([7]) for the kink profile and Eq. (l28l) for its tail 
(both with = 1 ), we find 


Xn =I2-\/2 [exp (—\/27„) — exp (—\/27„_i)] , [NEW approach] (39) 

in = 6-\/2 [exp (—V^in) — exp (—\/27„_i)] . [KO approach] (40) 

Equation (1551) agrees with Ei and Ohta m and with Carr and Pego [T^ . These authors use a perturbative approach 
where the small parameter is the extension of the domain wall defining the kink, but while Carr and Pego r^ on the 
existence of a Lyapunov functional, Ei and Ohta do not. Instead, Eq. (|40l) agrees with Kawasaki and Ohta Q, whose 
approach has been exemplified in Sec. Iml In Fig. [5] we compare old (dashed line) and new (full line) approach with 
exact kink dynamics (squares), showing that the new approach is quantitatively correct. 

• For TDGL4 equation, the two approaches give substantially different results, as we are going to show. In Eq. (1551) 
we now have only the term i = 2, with 04 = — 1 , and 


Xn = 


Er da: {dxKY 



+ U'\hJ 




(41) 


Now, see Eq. (1151) . R{£) = Acos{k£ + a) exp(—«;£), so that (even up to a constant) 


R{£) Y (0 and R{£) Y R"^ (0 • 


(42) 
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If we use the correct expression for R{tj we obtain 
2U" 

in = —r^ - - - [cos {Kin + 2a) exp {-Ojin) - COS {nin-l + 2a) exp {-Kin-l)] ■ (43) 

J_oo ^2^ 

In Fig.m we compare the full numerical solution of the continuum TDGL4 model (squares) with our results (Eq. (I43L 
full line) and with results obtained with the multikink approximation (dashed line). Our new approach of kink 
dynamics reproduces quantitatively very well the full numerical solution. In addition, the results from the multikink 
ansatz approach cannot be corrected using a simple rescaling of time, as in the case of TDGL. 


B. Conserved case 


Similarly to the nonconserved case, we are going to consider the general model 


dth = -dxx a2idfh - U'{h^ , 


(44) 


which requires more involved mathematics, whose details are partly given in Appendix Here we provide the final 
result. 


1 


'^h'^inin-l — An{in + in-l) 


{{in — 1 [in+lAn+1 /(^n+ 15 -^n — 1 )] 4” in [in —1 An—I f{in^in—2)]\ (45) 


where 


and 


An — 


n+4 


dx dxh / dx' {h — h^_i) 


= -flC'd)' 

-4 (I) Ri2^-2k) _ ^(2fc) 2[/" (I) - (I)) 




which reduces to 


4hlinin-l - 2V2{in + in-l) 


^in -1 2V2in+l + 8t/" (^exp{-V2in+i) - exp{-V2in-i)j 


d-in 


2V2xn-i + 8t/" (^exp{-V2in) - exp(--\/2Z„_2)) j | [CH] 


for the CH equation, and to 


I 


[CH4] 


dhiinin-l - A{in + in-l) 

[in-l [in+iA + 217" (cos {Kin+1 + 2a) exp {-nin+i) - cos {nin-i + 2a) exp (-k£„_i))] 
+in [in-lA + 2.U!/nA^ (cOS {nin + 2a) exp {-Kin) - cos {Kin-2 + 2a) exp {-Kin-2))] } , 


(46) 


(47) 


(48) 


(49) 


for the CH4 equation, with A = dx{h‘^ — h^). 

The previous two equations are rather involved and the expressions for kink speeds in are coupled, see the terms 
proportional to in±i on the Right Hand Side. Since the terms proportional to in in the Right Hand Side of Eq. (gSl) 
are smaller than the term in on the Left Hand Side by a factor ~ 1/in for large we may neglect them when 
in ^ a. Analogously, at denominators we can neglect the terms linear in i with respect the terms quadratic in i. 
Einally, we obtain a simplified version of Eqs. (I48I49I) : 


Xn — 


[/" (^exp{-V2in+i) - exp{-V2in-i)^ 

[CH simplified] 


dhl^inin-l 1 
+8inU// ^exp{-V2in) - exp{-V2in-2)^ 


(50) 
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and 

^ Ah'^ee —I ^ simplified] 

{in-i [2t/" (cos {nin+i + 2a) exp (-nin+i) - cos + 2a) exp (-k£„_i))] 

+£n \2U'^A^ (cos (k£„ + 2a) exp (-k£„) - cos {Kin -2 + 2a) exp {-Kin- 2 ))] } ■ 

In Figure [SJ we compare the different approaches and the numerical solution of the continuum CH equation, while 
in Fig. [6] we do the same for the CH4 equation. In both cases, exact numerical results are given by squares, our 
full analytical expressions Eqs. (148149^ are given by solid lines, our simplified expressions Eqs. (150151|) are given by 
dotted lines, and the analytical expressions using multikink approximations are given by dashed lines. The two figures 
clearly show that our full expressions (I48I49I) reproduce correctly numerics of the continuum model in both cases. 
The simplified model provides a reasonable result, but it is quantitatively inaccurate, proving that the subdominant 
terms ^ l/£„ are relevant for the interkink distances in used in the simulations of Fig. [S] and Fig. [51 However, these 
subdominant terms should become negligible for larger inter kink distances in- 


V. STABILITY OF STEADY STATES 

In Sec. imi we have shown that TDGL-kinks feel an attractive interaction while TDGL4-kinks feel an oscillating 
interaction, even if in both cases R{x) vanishes exponentially at large x. This fact implies two important differences: 
(i) all TDGL steady configurations are uniform, Xn+i — Xn = i, while TDGL4 ones may be even disordered; (ii) all 
TDGL steady states are linearly unstable, while TDGL4 steady states may be stable or unstable. Let us prove these 
statements. 

We can rewrite Eq. (1^^ incorporating the positive prefactor at RHS in t, 

Xn — R{Xn Xn—l) R{Xn-irl Xn)i (^2) 

whose time independent solution is R{in) = R{in-i) Vn, i.e. R{in) = x, with £„ = Xn+i — Xn- Eor the standard 
TDGL model R{x) is a monotonous function, so the equation R{in) = x has at most one solution. In practice, every 
uniform configuration £„ = £ is stationary. Instead, for the TDGL4 model, the equation 

R{in) = Acos{Kin + «) exp(—= r (53) 

has a number of solutions which increases when decreasing jrj, up to an infinite number of solutions for r = 0. 

As for the stability of a steady state, let us first focus on uniform configurations, i.e. all £„ = i. In order to study 
the linear stability of this configuration we need to perturb it, 

in{t) = i +en{t), (54) 

and determine the temporal evolution of the perturbations e„(t) <C i- Using Eq. (15^ we get 

en = 2R{in)-R{in+l)-R{in-l) 

= R'{£)(2en — e„_|_i — e„_i), 

whose single harmonic solution is en{t) = exp(wt + iqn), with 

w(g) = 4i?'(^) sin^ (I) . (57) 

We have stability (instability) if R'{i) < 0(> 0). Since R 2 {x) is an increasing function, see Eq. (^5)) . any uniform 
configuration is unstable for TDGL. This result leads to a perpetual coarsening dynamics Q. Instead, since Ri{x) is 
oscillating also its derivative is oscillating and with varying i we obtain stable steady states if R'{i) < 0 and unstable 
steady states if R'{i) > 0. 

In the general case of a nonuniform steady state, 

in{t) = in + in{t) with R{i*n) = r, (58) 

Eq. (1551) . which is still valid, gives 

= 2R\il)en - i?'(C-i)en-i - i?'(f:+i)e„+i. 


(55) 

(56) 


(59) 
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The linear character of the equations allows to write en{t) = getting 


2R'{OAn - - R'{t^+Mn+1 = aA 


(60) 


but the n-dependence of R'{(n) prevents the diagonalization with Fourier modes (A„ ^ 

Multiplying Ea. dBO)) with , summing aver all n, and after some simple recombinations of the l.h.s., we 

obtain 

N N 

^ (61) 
n—1 n—1 

which shows that eigenvalues a are real. Furthermore, if all quantities R'{^n) have the same sign, a has the sign 
of i?'(^*). In particular, any steady-state kink configuration with i?'(^*) < 0 for all n is stable. As a consequence, 
i?'(t'*) < 0 for all n is a sufficient condition for stability, and there is an infinite number of stable configurations in 
which the system can be trapped and stuck during the dynamics. 

If the quantities R'{(-n) exhibit both positive and negative signs, Eq. (ICT) does not allow to draw conclusions. 
However, in the simple cases of a period-2 configuration, = 1^+2^ ^ period-3 configuration, = £^+ 3 ^ '"^e can 

prove that R'{£n) < 0 is also a necessary condition for stability. Let’s show it explicitly for the period-2 configuration. 
If 


we obtain two coupled equations which are solved assuming 

A2„ = A2„+i = Cieh2-+1)9. 

The resulting eigenvalue equation is 

^2 _ 2 a + 4i?'(C)i?'(C) sin^ Q = 0- 

We have stability if both eigenvalues are negative, i.e. 

stability < 0 and R'{£* 2 ) < 0. 


(62) 

(63) 

(64) 

(65) 


VI. SUMMARY AND DISCUSSION 


Our paper studies kink dynamics deriving from a generalized Ginzburg-Landau free energy, see Eq. (1151) . The 
potential part of the free energy, U{h), is the classical, symmetric double well potential, typical of a bistable system. 
The “kinetic” part of the free energy is the sum of squares of order parameter derivatives of general order. 

The main motivation to study such free energy is that there are systems whose “kinetic” free energy is not given 
by surface tension, proportional to (h^), but rather to bending energy, which is proportional to (h^^.). Since the two 
terms are not mutually exclusive, it is quite reasonable to consider the free energy 


R = 


dx 




( 66 ) 


Then, we have generalized previous expression to Eq. m- However, even if our treatment is valid in full generality, 
we have focused on two cases: iFi = l,Ar 2 = 0 and Ki = 0 ,K 2 = 1, i.e. to pure surface tension systems (to check 
existing results) and to pure bending systems (novel system of specifical biophysical interest 0)- 

Once J- is given, we may derive a generalized Ginzburg-Landau equation, see Eq. m, or a generalized Gahn- 
Hilliard equation, see Eq. (|4^ . The standard approach to derive an effective kink dynamics is to assume a specific 
form of h{x,t) as a suitable superposition of kinks, hi,(x — Xn{t)), located in Xn{t). This method has proved to 
be fruitful, because it has allowed to explain coarsening dynamics of TDGL/CH models 0 [l3l - [T^ . to determine 
coarsening exponents, to study the effect of a symmetry breaking term [a , and the effect of thermal noise. 

However, the ability of the multikink approximation to reproduce quantitatively the exact dynamics of the contin¬ 
uum model was already questioned by Ei and Ohta m for the TDGL model. The failure of this goal is even more 
transparent when considering the bending energy, i.e. the TDGL4/CH4 models. In Figures [3][6] we make a detailed 
comparison of exact results (squares, derived from the direct integration of the equation) with the standard multikink 
approximation (dashed lines) and with our new results (full lines). The conclusion is that the new approach gives a 
reliable, discrete description of the exact, continuous dynamics: see how full lines follow squares in all Figs. [Sjni 
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We can still ask why we should derive an approximate kink dynamics if we have the full exact dynamics of order 
parameter h{x^ t). There are several good reasons: (i) an analytical approach to nonlinear full dynamics is hard if not 
impossible; (ii) kink dynamics is easy to understand and analytical methods are feasible; (hi) numerical simulation of 
kink dynamics is far faster than the simulation of the full PDE. 

In addition to be numerically reliable, some of our kink models (TDGL4/CH4) have the advantage of showing an 
oscillating tail R{x) = hi,(x) — h^. This oscillation implies two important features. Firstly, an oscillating tail means 
an oscillating force between kinks, as opposed to the classical TDGL/CH models. Therefore, the long term dynamical 
scenario is not a coarsening scenario, but the freezing in one of the many stable states Q- This can give rise to a 
consistency problem when we use the approximation ^ a to derive kink dynamics. However, the approximation is 
expected to give reasonable results even for not so far kinks and comparison with exact numerics supports such claim. 

Secondly, an oscillating tail R{x) is at the origin of a quantitative discrepancy between classical multikink approaches 
and our approach. Using numerical simulations, we have shown that our approach provides much better quantitative 
results. For example, classical results for TDGL4 provide an interkink force proportional to R{£) while a force 
F{£) Ri R?(ll2') appears to be more appropriate. If it were R{t) ~ exp(—/c£), the two approaches would be equivalent, 
apart a rescaling of time. Instead, if R{t) ~ cos(k.^ + a) exp(—the two approaches are definitely different. 

In this paper we have focused on the derivation of kink dynamics and on the quantitative comparison with the 
exact dynamics of the PDF. The kink models for TDGL4 and CH4 are also considered in Ref. [l^ where we specially 
use them for long time dynamics of the deterministic model and for any time dynamics of the stochastic models. In 
fact, once we have proven (here) their quantitative reliability, we can use them with confidence whenever the direct 
numerical integration of PDFs would be too demanding in terms of CPU time. This is certainly the case if we require 
to go to very long times or if we need to add stochastic noise to the equations. Our evaluation of the simulation 
times for the PDE (^pde) and for the kink model {F) allows to conclude that we gain four orders of magnitude, 

tpDE/tk ~ 10"^. 
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Appendix A: Derivation of Eq. (I45D 


Let us rewrite Eq. (HI) . 

dth = -a,, a2^^l^h - U'(h)^ (Al) 

and we still suppose we can write dth —Xndxh. If we integer (ED twice between oc^_i and x we obtain 

6.x' (h-/i„_i) = -^a 2 *a^*h + C/'(/i)+j„_i(a:-x„_i) + /r„_i, 

■-i i 

with j = dx (X^i a2id1^h — U'{h)) and fi = o, 2 idx^h — U'{h). Then we multiply by dxh and we integer between 
and : 


— 


k<^ 


^ dfhdf-^'^h I + u{h) 




n+T 


+ + Mn-iWnb 


(A2) 


where A„ = /Ei ^^h dx' (h — h^_i) and = /Ei dxh{x — x„_i). 
If we do the same thing but between x and x„_|_ i we obtain : 

(-l)-i 


^n-^n — 


k< 7 


- E ^ dfhdl^-^’^h + U{h) 




n+7 


n+i 


+ (^3) 
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where B'^ = -f^_i dx d^h(x„_^i — x). By summing Eq. 
kink we find 


for the n—th kink and Eq. (IA3I1 for the (n — 1)—th 




fc< A 


- E I ^ dfhdf-^^h I + U{h) 


fc=l 


n+^ 


(A4) 


+ + B^-l) + 

nn+i 


Because of the defintion of j, in[h]^_l = j„+i — jn-^- Finding j„_i from (|A4p and jn+^ from the same equation 
with n —>■ n + 1, we finally obtain the kink dynamics 

= 


Bn+1 + B'^ 

( 

XnAyi Xji+l-^n+l A 


V 


E I _ U{h) 


-I n+f 




- An+iW”!| 


1 


+ ^n-1 


/ 


n-l-i \ 

^n—l^n—1 ^n-^n H" 

^-^{dlhf - ^ dfhdl^-^^h - U{h) 


V 

_ i \ ^ k=l / 

^-1 7 


(A5) 


We now estabish the relations 


MT-1 = (-l)”2h„ 

|•n+^ /■"+5 

i?„+i + = / dx 9a;h(a; - - / dx 9a:h(a ;„+1 - a;) 

d„+i ^ Jn-^ ^ 

-1 


p+§ /•n+^ 

/ dx h„+3 + 

/n+i 


dx h^_i - 


dx h 


~ (-l)”+^2h„(a;„+i - x„). 

Then, reminding that = T.^a2i^fh |„+i — C/'(h„_i_i) is very small for a membrane near the equilibrium and 

[h] is very small in the limit of distant kinks, Eq. (IA5I) becomes 


-1 


4/l2 I Xn+l - Xn 


( 


^n-d^n d^n+l^n+1 A 


V 


E 


a2i 


(- 1 ) 


2-1 


k< ^ 


fc=l 


\ 

n+f \ 

- Uih) 

"-7 


^n—1 


/ 



in —1-^71 —1 in-^n 

^ a2. ^-^^{dlhf - ^ dl’^hdf-^'^h - U{h) 


V 

_ * \ ^ fc=i / 



(A6) 


or, collecting Xn terms in the right-hand side 

1 


Xn = 


{-^n —1 [Xn+l^dn+l A /(^n-t-1;-^n —1)] A in [^n —A f {in-; in— 2 )]} (14:5D 


dhlinin-l - An{in + in-l) 

whereAn and f{x, y) are defined in the main text respectively Eqs. (H51) and (ITTl) . 
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Appendix B: Numerical schemes 

We integrate the full dynamics of various PDE’s on a one-dimensional lattice with periodic boundary conditions. 
Space-derivatives are calculated using a finite-size difference scheme with discretization dx = 0.2. The time integration 
is performed using an explicit Euler scheme. The time-step depends on the equation to be solved. The most stringent 
case is CH6, where must be dt = 10“®. Therefore, we have used this value of dt for all simulations. 

Initial conditions are built from stationary kink profiles, which are known analytically for TDGL and CH. Eor the 
fourth order case, TDGL4 and GH4, these kink profiles are obtained numerically from steady-state solutions with 
isolated kinks. To build a kink-antikink pair at xi and X 2 , we use the profile of the stationary kink hi.(x — xi) for 
0 < a; < (xi -I- X2)I‘2‘, and the profile of the stationary antikink, —hi,(x — X 2 ), otherwise. 

For the implementation of kink dynamics we also use an explicit Euler scheme with dt = 10“^. The main difficulty 
comes from the annihilation of a kink-antikink pairs. Indeed, since kink models are designed to by quantitatively 
accurate at long inter-kink distances only, they are not necessarily accurate, or even well defined at short distances. 
However, kinks and antikinks merge and annihilate rapidly in the full dynamics when their separation is smaller than 
the size of the kink core. We therefore generically use a cutoff inter-kink distance a below which kinks and anti-kinks 
spontaneously annihilate in kink dynamics, using the following procedure. If the separation between the n-th and 
(n -I- l)-th kinks is smaller than a at time t + dt, we define the collision time-step dt^^^, = ixn+i{t) — x„{t))/{xn(t) — 
Xn+i{t)) from a simple linear extrapolation. We then integrate the dynamics of all kinks up to < -f dtanm and erase 
the two kinks at and Xn+i- 

In the non-conserved case, the full dynamics is not much affected by the choice of the cutoff, and we chose a = 0 
in the simulations presented in the main text. In the conserved case, the denominator (4:hm)‘^inin-i — + ^n-i) 

in Ea. (l4^ vanishes for positive £n and £„-!■ Assuming that £„ £n-i just before the collision, we find that the 
critical value of the interkink distance to keep the dynamics well defined is £c = 2/{2hmf' ■ Using the quartic potential 
U = —h^h'^/2 -I- h'^/4:, with hm = 0.9, we find A ^ 1.86, leading to £c ~ 1.15. We find that our numerical scheme 
is stable for a > 1.3, which is consistent with the expected constraint a > £c- In order to ensure strong stability, we 
have performed most simulations with a slightly larger value of a = 1.5. 
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